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Abstract: We extend the work of [1], including an imaginary chemical potential for quark 
number into the Sakai-Sugimoto model and codimension k defect theories. The phase 
diagram of these models are a function of three parameters, the temperature, chemical 
potential and the asymptotic separation of the flavour branes, related to a mass for the 
quarks in the boundary theories. We compute the phase diagrams and the pressure due 
to the flavours of the theories as a function of these parameters and show that there are 
Roberge- Weiss transitions in the high temperature phases, chiral symmetry restored for 
the Sakai-Sugimoto model and deconfined for the defect models, while at low temperatures 
there are no Roberge- Weiss transitions. In all the models we consider the transitions 
between low and high temperature phases are first order, hence the points where they 
meet the Roberge- Weiss lines are triple points. The pressure for the defect theories scales 
in the way we expect from dimensional analysis while the Sakai-Sugimoto model exhibits 
unusual scaling. We show that the models we consider are all analytic in when /x^ is 
small. 
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1. Introduction 



Since the original discovery by Maldacena of a duality between type IIB supergravity on 
AdSs X and = 4 super Yang-Mills at strong coupling [2] there has been a consistent 
effort to find a holographic model that is more like the theories we observe in nature, 
i.e. contain less supersymmetry and the correct degrees of freedom. In particular, much 
research has been dedicated to finding QCD like models at strong coupling using a classical 
or semiclassical supergravity. One such model that has seen some success is that of Sakai 
and Sugimoto [3] which is a holographic model in type IIA supergravity, dual to a QCD 
like large theory, where A^ is the number of colours. The boundary geometry is x 
so at suitable energy scales the gauge theory is effectively 4 dimensional, furthermore 
applying antiperiodic boundary conditions for fermions on the breaks supersymmetry 
so the only remaining degrees of freedom are those of large A QCD. Since its discovery, 
much progress has been made in understanding how quark number [4] and isospin [5, 6] 
chemical potentials affect the phases of the model, and also the effect of chemical potentials 
in holographic models in general [7, 8]. 

Another approach to studying field theories is to analyse them in different dimensions. 
It has been known for many years that the gauge theories in different dimensionality to 
the usual 3 + 1 we are familiar with can behave in very different ways than expected. 
Furthermore, much effort has gone into studying field theories in 1 + 1 dimensions, as 
a much more general class of theories can be solved exactly than in 4 dimensions. One 
may produce gauge theories in AdS/CFT with different dimensionality by altering the 
dimensionality of the branes that go into producing the background geometry. In particular, 
we will be interested in systems where the colour sector lives in 3 + 1 dimensions but the 
flavour sector lives on a 3 — A; + 1 dimensional submanifold of the boundary, dual to a 
background of A D3 branes with Nj D(7 — 2k) branes in the probe approximation added. 
The intersection of these branes forms a defect. It was found in [9] that the order of the 
confinement /deconfinement phase transition in the T/ fi plane depends on the codimension 
of the defect for real chemical potentials. 

The motivation for studying theories with an imaginary chemical potential comes from 
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the sign problem in lattice gauge theory, which makes simulating QCD with a real chemical 
potential hard to accomplish because the fermion determinant for such a system is complex. 
A way to circumvent this problem using standard lattice techniques is to take the chemical 
potential to be a purely imaginary number, and systems where this is the case have been 
extensively studied in lattice gauge theory [10]. 

The goal of this paper is to compute the phase diagram for the Sakai-Sugimoto model 
and for codimension k = 1 and k = 2 defect theories in the presence of an imaginary 
chemical potential. In section 2 we provide a short review of the analysis of Roberge 
and Weiss [11] who studied a field theory system with an imaginary chemical potential 
in perturbation theory. In section 3 we describe the results of [1] where the D3 - D7 
system, dual to = 4 supersymmetric Yang-Mills with N = 2 fundamental matter at 
strong coupling with an imaginary chemical potential was studied. Sections 4 and 5 give 
a short review of and the Sakai-Sugimoto model at low and high temperatures [3, 12, 13] 
and codimension k defect theories respectively followed by the computation of their phase 
diagrams. We conclude with a brief discussion and future outlook. 

2. Roberge- Weiss Transitions 

In 1986 Roberge and Weiss studied QCD with an imaginary chemical potential for quark 
number [11]. It is well understood that the phases of a pure SU{N) gauge theory at 
finite temperature can be distinguished by the expectation value of the Polyakov loop {P). 
At low temperatures the theory is in the confined phase which implies the expectation 
value of the phase of the Polyakov loop vanishes. This further implies there is a Zjv 
symmetry where the Polyakov loop transforms as P — )■ P with r an integer. For high 
temperatures this symmetry is spontaneously broken and hence {P) 7^ so the theory 
becomes deconfined. When one introduces fundamental matter into the theory the Z^r 
symmetry is never a symmetry of the theory, however a remnant of the Zjv symmetry is 
still present. Introducing an imaginary chemical potential makes this explicit and resolves 
the ambiguity between confined and deconfined phases of the theory. Consider the partition 
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function for such a theory 

Z = TV ('e-/^-f^+^/^w^«) (2.1) 



Nq is the quark number operator and /i/ is the imaginary chemical potential^. One can 
see this partition function has a symmetry under the transformation |U/ — )• yUj + ^ where 
/3 is the inverse temperature of the system. Together with the fact that A''^ is quantised, 
this imphes G [— vr, vr]. 

The adjoint sector also contains a Z^v symmetry which introduces a further periodicity 
into the theory. The imaginary chemical potential couples to fundamental matter in the 
same way as the time component of the SU{N) gauge field. Hence, the time component 
of the gauge covariant derivative associated to this gauge field is modified by the presence 
of this chemical potential in the following way 

dr + iAr dr + lAr - ifij (2.2) 

We can remove the dependence of the action on ^u/ by rescaling the fields so that the 
dependence on /i/ is in the boundary conditions around the thermal circle 

(/) (x, ^) ~ (x, 0) (2.3) 

where the fermionic and bosonic modes differ by a sign. We may now apply an SU{N) 
gauge transformation with the group element U (x, r) with the property that U (x, j3) = 
e N U (^x, 0) and r an integer. The action and path integral measure are left invariant but 
the boundary conditions around the thermal circle are not, hence 

^{x,P) ^ e<'"'^+'-^)(j){x,0) (2.4) 

The partition function must be invariant under these transformations, therefore 

27rr 



Z [fii] = Z 



(2.5) 



which further implies I3^i G [~;^; • 

The free energy F [fij] = — Tlog[Z] of QCD in perturbation theory, valid at high 
temperature, was first computed in [11] who found first order phase transitions at constant 
^Note that we have made the replacement ii — >■ ijii so fj,i on it's own is in fact real. 
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values of fij . These values are 



(2r- l)7rr 



(2.6) 



Lattice studies suggest that at low temperature when the theory is strongly coupled there 
are no phase transitions and the physics is smooth as a function of /i/, while at high 
temperature there are first order Roberge- Weiss transitions at m = ^urw [10]- 

3. The D3-D7 System 

It is possible to study A/" = 4 super Yang-Mills with M = 2 fundamental matter at strong 
coupling with an imaginary chemical potential analytically via the celebrated AdS/CFT 
correspondence. In order to introduce a chemical potential into the D3 - D7 model it is 
necessary to introduce a non zero background time component of a worldvolume gauge field 
on the D7 branes. Furthermore in the Euclidean formulation to get an imaginary chemical 
potential this constant value of the gauge field must be real. The quark number chemical 
potential fi of the boundary gauge theory is given holographically by 



At{r) is the t component of a U{1) gauge field and r is a radial coordinate where r — )• oo 
is the boundary limit. The gauge field strength is F = Ffj^^dx^'' A dx'^ and since we do not 
wish to study the theory as a function of a conductivity the only non zero component is 
Ffrdr A dt. Since we work in Euclidean signature one can see that the gauge field must be 
completely real in order that the the chemical potential is purely imaginary since when we 
Wick rotate the time direction the field strength becomes iFrr-dr A dr. The features of 
the phase diagram are as follows; At low temperatures there are no phase transitions for 
varying chemical potential. The first order meson melting transition that appears for zero 
chemical potential extends into a first order line at T = Tc (/x/). The dependence of Tc on 
/i/ was found to be 



Where Tq = Tc(0) ~ 0.77, K ~ 33.5, A is the 't Hooft coupling and m is the quark 
mass. Q is the expectation value of the phase of the Polyakov loop. Higher corrections 



^ = lim At{r) 
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Figure 1: Phase diagram of the D3 - D7 system with an imaginary chemical potential. 



are suppressed by factors of 4 since 



Tc{0) 



a 



< Jr. Above this temperature there is 



an infinite set of Roberge- Weiss transitions at constant values of /x/ = /xrw = 

These Roberge- Weiss lines are first order and so the points where they meet the meson 

melting line are triple points. 

The argument of [14] (described in a later section) leading to the presence of Roberge- 
Weiss transitions in gravity duals is quite general. We expect a similar structure in both 
the defect theories we consider and the Sakai-Sugimoto model from the argument given in 
[8]. The details of the theory such as the Tq = Tc(/i/ = 0) and the phase transition that 
takes the place of the meson melting transition here will not be the same however. 



4. The Sakai-Sugimoto Model 

In order to model QCD at strong coupling holographically the Sakai-Sugimoto model takes 
a stack of N D4 colour branes in type IIA supergravity and in the large N limit replaces 
them with a background geometry. To incorporate Nf flavours D8 - D8 branes are added 

^One might imagine that, since the spacing of the Roberge- Weiss phases is of order l/N, we should 
consider a model that goes beyond the supergravity limit, however, tlie factors of arise from the symmetry 
of the theory under shifts of /ij — + so it is consistent to consider the supergravity limit only. 



-6- 



which, if Nf <^ N we can consider in the probe approximation. To reduce the dimension- 
ality of the boundary from 5 to 4 one of the dimensions is compactified on a circle, so if 
we consider meson masses below the Kaluza Klein mass scale Mkk the boundary gauge 
theory is effectively 4 dimensional. In addition, imposing antiperiodic boundary conditions 
for fermions on the compact direction breaks supersymmetry completely making the only 
light degrees of freedom those of large N QCD. A great successes of the Sakai-Sugimoto 
model is it is a holographic model that realises chiral symmetry breaking. The geometry 
resulting from the D4 branes contains a horizon, and the D8 - D8 solutions that end at 
some finite height above that horizon represent a state in the field theory where chiral sym- 
metry is broken. If one increases the temperature with the other parameters of the theory 
fixed the distance from the horizon to the end point of the D8 - D8 branes is reduced. 
Eventually, the horizon will meet the branes and chiral symmetry is restored. It has been 
shown in [12] that there are two bulk geometries with similar asymptotics, corresponding 
to QCD at low and high temperatures. For low temperature 



3 , , 3 



ds^ = (I) ' {dtl + 6,,dx'dx' + fKK{U)dxl) + (^1) ' {j~(U) + ^''^^') (^-^^ 



with fKK{U) = 1 — i^^W^^ ■ R is the radius of curvature and is related to the string 



coupling and string length. C/kk is the position of the horizon and both are related to the 



Kaluza Klein mass Mkk = which is the mass that characterises the scale of the 

compact direction. For the high temperature deconfined phase we have 

3 3 

ds^ = (I) ' {fT{U)dtl + 5,,dx'dx^ + dxl) + ) ' + U^dnl^ (4.2) 

where friU) = 1 - (^ff-^ . Note we have Wick rotated to Euclidean signature, and to 
avoid a conical singularity at the horizon the Euclidean time direction must be periodically 
identified, hence there are two compact directions in the geometries and thus one can think 
of the low / high temperature behaviour as a competition between which circle shrinks to 
zero size first. Ensuring regularity at the horizon one can deduce that the temperature is 
given by 
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We will consider both background geometries, but we will pay particular attention to the 
high temperature phase because we expect to see Roberge- Weiss transitions there. In this 
case the tE circle shrinks to zero size, and thus the spacetime ends, at U = Ut while the 
2:4 circle remains finite sized there. 

The phase diagram of the Sakai-Sugimoto model will be a function of 3 variables, 
the temperature and chemical potential are obvious but there is also another, less evident 
parameter. In the gravity picture there is a parameter L that is the asymptotic separation 
of the D8 - D8 branes in the X4 direction. Since in the effective field theory this direction 
is compactified and integrated out, operators like ^ip^^ A^iptp^'^ Ayij: will be generated with 
coefficients that depend on L. These operators are irrelevant, so it is not clear exactly 
how to relate L to the various generated operators when the coupling is large. We will 
show below that the temperature and chemical potential naturally combine with L to form 
dimensionless parameters, which we will use to plot the phase diagram. 

4.1 D8 Embedding 

Consider the metric for the D8 branes in the background (4.2). We choose the brane to 
share all coordinates except the direction where the D8 is pointlike, which we make a 
function of the radial direction U . The induced metric of the D8 with this choice is 

dsl, = (^'" {fTiU)dtl + 5,,dx^dx=) + [{^^y + (I) ' ^4 ) dU'+ 

3 

+ (^^yU^dnl (4.4) 

Prime denotes differentiation with respect to U . To incorporate a chemical potential into 
the gauge theory, on the gravity side we include a worldvolume gauge field on the D8 
branes. We also introduce a NS B field into the bulk, dual to the expectation value of the 
phase of the Polyakov loop in the gauge theory. The only non-zero components of these 
field strengths are Ftu and Btu- We can now write the DBI action of Nf branes 

5dbi = iV/TD8 j d^xe"-^ Vdet {gab + Bah + 2TTa'Fab) (4.5) 
= / rf^xdC/C/y ^ (1 + iB,u + 2na'F,uf) + fTiU)x'i (4.6) 
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Td8 is the brane tension and we have used that the dilaton is given by 

Furthermore, we may set the U components of the gauge potentials to zero with an appro- 
priate choice of gauge, therefore Ftjj = —dijAt = —A' hence Chern Simons terms do not 
contribute to the action. 

Since the action only depends on derivatives of X4 and A we have very simple equations 
of motion that after one integration take the following form; The equation of motion for 
Xi is given by 

'fT{U)x'i + {§f (1 + (Btu + 27ra'Ftuy 

and the equation of motion for Ftjj is 

UR^ {Btu + 2Tia'Ftu) dg^ 



. ^ ,0 /Rn3/ . .2\ 2TOWfrD8Vol (54) 



(4.9) 



where k and d are constants determined by boundary conditions. There are two different 
types of brane profile allowed by these equations of motion. The first is where the brane 
extends all the way to the thermal horizon, covering all of the radial direction, and secondly, 
where the brane turns around meeting a D8 at some finite radius Uq. In the low temperature 
phase only the latter solution is allowed, because the X4 circle in the low temperature 
geometry shrinks to zero size first while the thermal circle remains finite sized over the 
whole range of the radial coordinate. The converse is true in the high temperature case, 
allowing both types of solution described above. These solutions are dual to the chiral 
symmetric and broken phases in the boundary theory respectively. 

It is straightforward to fix the constants above by considering boundary conditions. 
Consider (4.8) at some point such that x'^ (U) = 0. We find 

k = (X4) (4.10) 

Therefore k is forced to vanish. We will find a useful constraint on d for this type of 
embedding later. 
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Consider now the other type of solution where the brane turns around. There must be 
some point U = Uq where the slope of X4^(U) diverges so consider x'^{U = Uq) — t- oo. From 
(4.8) we find 

k = U^VfTm (4.11) 



and from equation (4.9) 



d = 0{^] (4.12) 



So for this type of embedding, d is forced to vanish'^. Note: This analysis is analogous to 
the real chemical potential case studied in [4]. 

Consider the case where the D8 branes extend all the way to the thermal horizon; from 
equation (4.9) we can solve for {Btu + 2TTa'Ftu) and plug back into (4.8). Doing this we 
find 



and hence 



(4.14) 



One can see that (4.14) admits constant solutions that cover the entire range of the U 
direction, falling into the horizon at U = Ut- Since for these solutions 2:4 is arbitrary 
and any value of x^ will be have the same action, these solutions must correspond to 
the phase where chiral symmetry is unbroken and the quarks in the boundary theory 
are exactly massless. To ensure the reality of {Btu + 2Tra'Ftu) to produce an imaginary 
chemical potential in the gauge theory, and to ensure the reality of X4 one can see there is 
a constraint on d from (4.13) 

< R^U^ (4.15) 

In contrast to the condition derived in [1], the condition here does not depend on the 
horizon value of X4{U). In the D3/D7 case the bound on d depended on the asymptotic 



^This is strictly d = O j which does not necessarily constrain d to vanish if Fut diverges. One may 
assume U — U{x4) rather than X4 — X4{U) as used to derive (4.4) to show that this is not the case and d 
does indeed vanish for this type of embedding. 
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value of the slipping mode, or the quark mass in the language of the gauge theory. As 
described above this type of solution only corresponds to strictly massless quarks. 

4.2 The Potential 

4.2.1 The unbroken chiral symmetry case 

Firstly, we note that for all phases of the Sakai-Sugimoto model the Lagrangian quark mass 
is zero, furthermore in this phase there can be no dynamically generated quark mass as that 
would break chiral symmetry. As mentioned previously it was shown in [1], the chemical 
potential is related to the worldvolume gauge field on the gravity side while the expectation 
value of the Polyakov loop in the theory is related to the NS-NS 2 form. Together we have 

a-/3M/ = ^^ {Btu + 27ra'Ftu) (4.16) 

where the integral is over the thermal cigar. Plugging into this using (4.13) we find an 
explicit expression for a — 

a-/3fii = -^ dU (4.17) 
27ra' Ju^ y^U^RS _ ^2 

To find the potential we evaluate the (renormalised) action and eliminate d to find S in 
terms a — f3fj,j . In the simple case of black hole embeddings the integrals can be done 
analytically 



d^ 

(a - 



2-Ka'T 



3 



(4.18) 

Plugging into the action (4.6) with the equations of motion and the constraints found for 
black hole embeddings above we find 



^ ^ VWNfT^,\o\ {S') Vol (M3) /3 /- 



9s JUr _ g 

Let us relabel the overall constant A = ^^^f'^°^^°K^ )Vo1(m )/3 ^ This integral is divergent, 
SO after regulating we may add a counterterm to cancel the divergence 

S = A dU , - A^— (4.20) 
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Evaluating the integral and taking the limit A — )• oo we find 



1 1 6 fi'C/j, 
' 5 ' J2 



dp 



(4.21) 



Using d as a parameter we find figure 2. Since the potential is truncated at order it 




Figure 2: The Robcrge- Weiss potential as a function of a — /3^/ 

is instructive to expand these expressions for small d and eliminate to find 5 in terms of 
a — we find the first order correction to the action as a function of (a — /3^/) 



a — fipii 



d 



+ 0[d 



3Tra'TR2U^ 



(4.22) 
(4.23) 



hence 



S D 3TT'^a'^AT^U}{a - /3fiif + O {{a - (3^)^) 



(4.24) 



We wish to write this expression in terms of field theory quantities; in the high temperature 
phase the X4 direction is finite sized over the whole range of the radial coordinate so there 
is no constraint on it to ensure regularity at the horizon. The size of the X4 direction is 
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therefore related to the Kaluza Klein mass scale simply by 

1 



27ri?,, (4.25) 



Mkk 

The tension of the D8 brane is a dimension 9 object related to the string length by Tds = 
(27r)® . We can relate the radius of cm'vature and the string coupling to gauge theory 
quantities using the following (from [12]) 

f <^-^^) 

Where we have used that 

gl = = MkkqI (4.27) 

Using the relation between the temperature and the position of the horizon (4.3) we find 
the potential due to flavour degrees of freedom in terms of field theory quantities is 

We have used that Vol (S"^) = j^. The potential from the flavours is quadratic, and we 
expect the potential due to the adjoint degrees of freedom to be periodic from an argument 
in [14] for type IIB, generalised to type IIA in [8]. When the theory contains only adjoint 
degrees of freedom there is an infinite set of degenerate minima in the potential as a 
function of a. This degeneracy is lifted when one includes fundamental matter and an 
imaginary chemical potential in the theory. The potential for the flavours only depends on 
the combination a — pfij so if we shift fij by ^ the shift can be reabsorbed by a — )• a + 27r, 
hence is only defined up to a factor of Schematically the effective potential is 

V,s = Va + Vf^ (mm(a-^) +^{a-Pfii)A (4.29) 

This is shown in figure 3 for only adjoint degrees of freedom and for adjoint and funda- 
mental degrees of freedom with zero chemical potential. For the latter it is clear there is 
a single global minimum and many local minima. As we increase the chemical potential a 
competition develops between the global minimum and an adjacent local minimum until 
they become degenerate, which happens when ^/ = where /iRw is defined in (2.6). 
Increasing the chemical potential further results in a first order transition from one minima 
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Figure 3: The effective potential due to the adjoint degrees of freedom (left) and both adjoint and 
fundamental degrees of freedom (right) as a function of a. 

to the next (see figure 4). In the new vacuum the expectation value of the phase of the 
Polyakov loop differs from the previous vacuum by 




Figure 4: The effective potential for adjoint and fundamental degrees of freedom when fij w /xrw 
(left) and fij > /xrw (right) 



4.2.2 The broken chiral symmetry case 

In the language of the gravity theory this corresponds to Minkowski brane embeddings 
that turn around at some finite Uq > Ut- In contrast to the previous case, there must be 
a dynamically generated quark mass in this phase, which is given by 

Hence, chiral symmetry is no longer a symmetry of the theory. Moreover, as shown 
previously d must vanish for these solutions. Together with (4.13) this implies that 
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Btu + 2-Ka'Ftu = and hence the potential does not depend on a — /3/i/. Since /x/ is 
interpreted as limij^oo {U)] we are still able to introduce a chemical potential via the 
introduction of a constant gauge field, but the physics remains smooth for all values of this 
chemical potential. 

This analysis is also applicable to the low temperature geometry where this type of 
D8 - D8 embedding is the only one allowed. The conclusion is the same and the potential 
does not depend on a — /3fij . 

4.3 The Phase Diagram 

At low temperatures the only solutions allowed are Minkowski embeddings so the phase 
diagram has no phase transitions. Above a certain temperature there is a transition from 
(4.1) to (4.2) and both types of solution are then allowed. We will compute the phase 
diagram of the theory at high temperature by considering the difference in the DBI actions 
of the two types of solution, dual to the grand canonical ensemble in the gauge theory. 
After applying the conditions for each of the solutions derived above and substituting in 
using the equations of motion we find 



'Mink 



/i^SiV/TDsVol (S^) Vol (m3) /3 /■°° ,,,,,5 U^f{U) 



hence we define 



= - (^Mink — 'S'bh) 



oo 



U^f{U) 1 \ /-^o u 



Uo 



^ ^ ^/fi^iV/^D8Vol {S') Vol (R3) p 

gs 

When AS" = there is a phase transition. It is worth making some remarks on the method 
used in [4] to compute the phase diagram of the Sakai-Sugimoto model in the presence of 
a real chemical potential. The following dimensionless coordinates were used to eliminate 
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Uq from AS 



ut 



Ut 



U 



(4.35) 



where Wick rotating c would recover the real chemical potential result. In terms of these 
coordinates we have 

' ~ 0,5 



AS- 



Dimensionless 



duu 



1 



u^fiu) - /(I) ^ 



du 



(4.36) 



Henceforth we shall refer to the difference in the actions in dimensionless coordinates of [4] 
as A5Dimcnsionicss5 while in the dimensionful coordinates we shall denote is AS. Evaluating 
(4.36) for c = or for zero chemical potential we find figure 5. There are two zeros, one 
AS 



0.02 



-0.02 



-0.04 - 




Figure 5: The plot of AS against ut — for zero chemical potential 

of which occurs at approximately ut = 0.74 which corresponds to the fi = limit of the 
real chemical potential phase transition line. This was originally computed in [12, 13]. The 
second is at ut = 1 which, when one considers a real chemical potential is lifted and does 
not contribute to the phase diagram. When we consider an imaginary chemical potential 
however, this zero persists and would lead to a second phase transition line on the phase 
diagram. We will show that the phase transition line that this zero leads to should not, in 
fact be present on the phase diagram. Consider the expression relating the dimensionless 
horizon position and the temperature (from equation (4.3)) 

' ^ (4.37) 



4tt V i?3 
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The asymptotic separation of the D8 - D8 pair is L and is given by; 



where 



L = 2 lim Xi{U) 

l-QO 

Juo 

/ dU 



[/8 f{U) 



fm 




F{ut) 



F (ut) = 2 



du- 



1 



(4.38) 



(4.39) 



Note that for black hole embeddings L is arbitrary, so we consider Minkowski embeddings 
by plugging in for x'^ using (4.8). Using (4.38) to eliminate Uq from (4.37) 



TL = (ut) 



(4.40) 



TL as a function of ut is not a monotonically increasing function and has a global maxi- 
mum, it is shown in figure 6. Suppose we consider a fixed T and L. There are two allowed 

TL 
0.15 



0.10 



0.05 



Uj 




Figure 6: TL as a function of ut 



Minkowski embeddings that have the same value of TL but have different values of ut- For 
a fixed ut however there is only one Minkowski embedding with a particular allowed value 
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of TL. We wish to compute the phase diagram by considering the action for all allowed 
solutions at each TL and chemical potential (there are three - one black hole embedding 
and the two Minkowski embeddings), the one that has the lowest action dominates. Using 
ut as a parameter to do this is problematic because for each ut we are only considering 
one of the Minkowski embeddings allowed at that particular TL. Computing the phase 
diagram using a parameter that only considers a subset of allowed solutions leads to a 
phase diagram with a phase transition line that would not appear were we to consider the 
full set of allowed solutions. 

The analysis of A5Dimcnsionicss has one important consequence however: we note that 
A'S'Dimensioniess and the Combination TL are functions of uj- only (at zero chemical poten- 
tial), therefore the action is a function of TL, rather than T and L individually. Consider 
now the chemical potential. Rewriting (4.17) using the dimensionless coordinates we have 



hence, the combination of the chemical potential and the expectation value of the phase 
of the Polyakov loop is a function of ut and c, and vanishes when c — )• 0. There are two 
scales on the gravity side (c and ut) that correspond to two meaningful scales on the gauge 
thoery side; TL and -^j^ (a — /3/i/), which we will denote by T and JL respectively. We will 
plot the phase diagram of the model as a function of these two parameters. The phase 
diagram is computed by finding the zeros of AS (4.33) as a function of Ut, d for a fixed C/q 
and converting these coordinates to T and JZ. We note that repeating this computation for 
various values Uq reproduces the same phase diagram, which is shown in 7. The blue line 
is the first order phase transition line that extends from the phase transition point found 
in [12, 13] for zero chemical potential. The shaded area of the plot is the area where the 
chiral symmetry broken phase dominates and the unshaded area is the one where chiral 
symmetry is restored. 




(4.41) 



eliminating Uq as before we have 




(4.42) 
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Figure 7: The phase diagram of the Sakai-Sugimoto model with an imaginary chemical potential 
as a function of the dimensionless temperature and chemical potential. 

The issue we have noted regarding the use of the dimensionless coordinates to study 
the phase structure of the model was not probelematic for real chemical potential. We 
will now show why this was the case. Consider A5Dimensioniess for the theory when Jl is 
small. Expanding in powers of c and similarly expanding equation (4.42) we may capture 
the leading behaviour of A^Dimcnsionicss as a function of Ji 

9jf 



AiSj^imcnsionless ~ AiSDimensionlcss 



A ^Dimensionless 



c=0 



c=0 



52F {utY \Iu\ 



52T 



(4.43) 



where we have used (4.40) to exchange F (ut) for T. AS'Dimcnsionicssl^^Q must be evaluated 
numerically and is shown in figure 5. The factors of ut in the denominator of the above 
expression are problematic at first sight, however, we note that the limited range of ut 
from (4.15) implies Uj< > (? so for a fixed value of c > 0, ut must be strictly greater than 
zero and the first correction term never diverges. When ut is close to 1, A5Dimcnsionicss|^^o 
is small and positive. For an imaginary chemical potential the first correction term that 
appears at order Ji^ is negative leading to the second zero of AS'Dimensionioss- Were we to 
Wick rotate Ji so that we are considering the theory with a real chemical potential the sign 
of this correction term would be positive, so ASoimcnsioniess does not vanish when ut is 
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close to 1 for a finite but small real JL. 

Performing a similar expansion for small d on AS (4.33) is also instructive. Using an 
expansion of (4.17) to eliminate d and (4.3) to eliminate Ut we find 

A5 f« A5| " y ~ ^^^^^ vr^iJsT^a'^ + O (^{a - /3/U/)^) (4.44) 

We could in principle restore the dimensionful coefficient and rewrite AS using field theory 
quantities, however it is not very enlightening to do so. It is important to note that the first 
correction term does not depend on Uq and the equations relating T to Ut and (a — /3/i/) to 
d contain no dependence on Uq so the above expression for AS contains no hidden factors 
of Uo- Therefore the system must be well behaved when Uq approaches Ut, which when one 
considers the model using the dimensionless coordinates discussed above corresponds to the 
Ut ^ 1 regime. This analysis also shows the action depends on the square of the chemical 
potential, ensuring the analyticity of the phase diagram as a function of (a — /3^/)^. 

In summary, when comparing the actions of the two different types of embedding it 
is important to work in a coordinate system where there is a one to one correspondence 
between the position of a horizon and a temperature. 

In the low temperature phase when T < where is the deconfinement temperature 
chiral symmetry is broken and the glue in the theory is confined. Above the deconfinement 
temperature in the chiral symmetry broken phase the glue is deconfined and the theory 
contains mesons. In the chiral symmetry restored phase the mesons melt and the fun- 
damental degrees of freedom in the theory are quarks and gluons. The quarks are also 
massless due to the chiral symmetry restoration. Including the Roberge- Weiss lines the 
phase diagram of the theory is given in figure 8. The equation of the phase transition line 
is 

T = 0.154 + 0.447^2 + (/I^) (4.45) 
The equation for T behaves as expected, ie, is a function of even powers of p. 

4.4 Pressure due to Flavours 

We will compute the pressure holographically by studying the DBI action as a function of 
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T 



Figure 8: The phase diagram of the Sakai-Sugimoto model with an imaginary chemical potential. 
The black line extends from the line found for the real chemical potential case. The red lines are 
the Roberge- Weiss lines which appear in the chiral symmetry restored phase. Since the equation 
of the phase transition line is a function of JI^ the second derivative of the phase transition line is 
strictly positive so even though the Roberge- Weiss lines appear at order the chiral symmetry 
breaking transition does show curvature between the Roberge- Weiss lines. 

temperature. From basic thermodynamics 



F = TdS - PdV + ndN 



(4.46) 



The pressure is given by 



P 



dF 
dV 



(4.47) 



holographically, the DBI action is dual to the free energy density in the grand canonical 
ensemble, therefore it is also related to the pressure. The DBI action is 



/ 



(4.48) 



The integrals over the field theory directions contribute a volume of and a factor of /3, 
because £dbi is not a function of these directions. One may go from the energy density to 
the pressure by dividing through by these factors or we may directly compute the pressure 
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by evaluating 

P = - j dUdn^CBBi (4.49) 

4.4.1 Black Hole Embeddings - Zero Density 

When the density is zero the pressure takes a simple form 



oo 

5 



Pbh = -^f / dUU2 (4.50) 



we have relabelled the overall constant in front of the pressure as 



^^^tVWm (4.51) 

Regulating the integral and adding a counterterm to cancel the divergence as A — )■ oo we 
have 



, 7 



Pbh = - lim TV 



2U2 



A— !>oo \ 7 

Hence 



^ 2Ai' 



(4.52) 



Pbh = yAAf/l (4.53) 



Using the expression for the temperature given in (4.3) we can show 

ftH^f<^?^fl"T' (4.54, 

7 3«5s 

Using the expressions given in section 4.2 we find the pressure in terms of field theory 
quantities is 

^™ = ^^^(mk^) <"^' 

where A = g\N . This expression is quite surprising, in particular the pressure is propor- 
tional r^. 

4.4.2 Non Zero Density 

Expanding Pbh for small d we have 

Pbh = M £ dU + + O (d^) \ (4.56) 



Evaluating the integral as before, showing only the first order correction for non zero d 

2?d? 



Pbh D -AA ( ^^^^ + O (rf") 1 (4.57) 
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For small density we can expand the expression for the chemical potential (assuming the 
phase of the Polyakov loop vanishes) 



1 Z""" d 

Hi = -- — 7/ dU , = (4.58) 
27ra' Ju^ ^/U^R^ - d? 

,,^S_^(?I^fM«aV„^0(d) (4.59) 



where we have restored the correct dimensionality of the density. In terms of field theory 
quantities, the density in terms of the chemical potential is 

The first order correction depends on the square of the chemical potential and it is 
Together, the leading order behaviour of the pressure is 

Since the pressure is a mass dimension 4 object, the expression for the pressure will always 
be proportional to multiplied by a dimensionless function of the scales in the problem, 
in this case T, fij and Mkk- An expected feature of the expression for the pressure is that 
the leading order term dominates when the temperature is very large but, as commented 
previously, the pressure computed here does not correlate with what is known in QCD. The 
action for the black hole embeddings can be evaluated exactly in terms of a hypergeometric 
function, and it is 

3 



P = 



7.22387ri(r-l)rH V^kk; V ^^^^ 

+ ^(T-l)HhFi (^,l,^,r) j (4.63) 



where 

^~ (27r)2 24312 M4j,d2 



^ ~ \2 o4q12 )i/r4 M (4.64) 
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Since this formalism captures the strong coupling behaviour of the field theory expanding 
for large 't Hooft coupling is a natural thing to do 

P^T^(^^j:L_J^^f^ + o(^]] (4 65) 

V 2vr 7.38 M3j^ iV^iV7A \xy J ^ ' 

When one substitutes in for small d in terms of m using (4.60) the same result as for the 
small density expansion is obtained. 

4.4.3 Minkowski Embeddings 

The integral for the action for Minkowski embeddings cannot be evaluated in closed form. 
The pressure for the Minkowski embeddings is 



'"''^ = -''C/""l u^nu)%m (4.66) 



Rewriting: 



Replacing Uq using 



W = -AA/ dU (4.67) 

V ^ uf{u) 



Ina'mg = M = Uq - Ut (4.68) 



and expanding in powers of M (note, since the lower limit of integration is a function of 
M we expand the integrand, evaluate the integral and then expand the result. The first 
order correction only receives contributions from the first order of the expansion of the 
integrand, while the corrections at order contain contributions from higher orders of 
the expansion of the integrand.) 

In the zero mass limit the pressure for Minkowski embeddings captures the same 
behaviour of the chiral symmetry restored phase so we choose to only consider the terms 
proportional to M as the leading behaviour. To O (M^) we have 

^iiink = -M \ d\J\ — —^^^5 ^ + ^ + O (M^) 

(4.69) 

Once we have performed the integral the leading behaviour of the pressure is 

Puink = -^MU^^^ I 5 log I 3.2^ J + 5V37T - 12 J (4.70) 
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using the above expressions to rewrite in terms of field theory quantities we have 



Mink 



A2 NfN 4 m„T / , / 64A2 \ \ ^ ^ 

T-s^—T^—^ I 5log\ ^ — ^ +5\/37r-12 (4.71 

5.2436 27T M^^[ ^[ 27 m^^M^^J j ^ ^ 

therefore the leading behaviour is proportional to J?^ log ( - 



4.4.4 Hadron Resonance Gas Model 

Since the results we have here are for a system which contains bound states of quarks and 
strong coupling it is natural to compare to the hadron resonance gas model of QCD [15]. 
This is a model of non interacting hadronic and mesonic resonances, which is a good model 
of QCD at low temperatures since the quarks are confined and are no longer the physical 
degrees of freedom of the theory. The pressure contribution of a single particle in this 
model is given by 

= I2 E (-^)'^' (fe/3m.) (4.72) 

k=l 

(3 is the inverse temperature, gi is a degeneracy factor and rj = 1 for bosons and —1 for 
fermions. K2 is a modified Bessel function of the second kind. Assuming the mass is small 
we have 

k=l ^ ' 

evaluating the sum 

fi^P, = ^ f-2Li4 (-r?) + 4^Li2 (-r?) + O [m\)\ (4.74) 



Note that Li is a polylogarithm function that contributes a numerical factor when we plug 
in for T]. Note that the first order correction to the pressure for non zero quark mass is 

P-T^% (4.75) 

While we note that the Sakai-Sugimoto is very different from QCD, being 5 dimensional 
and having a large number of colours etc., the pressure from the hadron resonance gas 
model exhibits very different behaviour to that found in the Sakai-Sugimoto model above. 
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5. Defect Theories 

The work of [1] showed that when one considers = 4 super Yang-Mills theory with 
quenched M = 2 flavours at strong coupling, dual to type IIB supergravity on AdSs x 
with probe D7 branes, in the presence of an imaginary chemical potential the resulting 
phase diagram contains only first order phase transition lines. Therefore, the points where 
the Roberge- Weiss lines meet the deconfinement transition are triple points. It would be 
interesting to find a model where the deconfinement line is not of the first order and one 
way that could possibly be the case is in a defect theory. This type of background was 
discussed in a very general way in [9] , but we limit ourselves to the background that comes 
from the backreaction of D3 branes in type IIB supergravity. To include fiavours we put 
Nf D(7 — 2k) branes into the AdSsX background, k is the codimension of the defect 
which we take as A; = 1 (A: = 2) representing a bulk theory containing Nj D5 (D3-D3) 
probe branes in the geometry that results from N D3 branes where — )• oo. This bulk 
theory corresponds to a field theory where the fundamental matter lives on a 2 + 1 (1 + 1) 
dimensional submanifold of the 3 + 1 dimensional boundary. The adjoint fields live in the 
full 3+1 dimensions of the boundary. A; = would correspond to having no defect, ie the 
D3 /D7 system. In order for the theory to be retain M = 2 supersymmetry, the brane must 
cover 2) — k spatial dimensions of the boundary and ?> — k directions of the in addition 
to the time direction and the radial direction. This type of theory has been discussed 
extensively in the literature. See for example [16]. 

The metric corresponding to AA = 4 super Yang-Mills at high temperature having Wick 
rotated to Euclidean signature is 




(5.1) 



with 




(5.2) 



Rescaling the radial coordinate using 



da 



dr 



(5.3) 



a 



r 
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leads to the metric 

2 



ds' = {^fh+{a) i dr' + 6^,dx^dxA + [da' + a'dQl) (5.4) 



where i and j run from 1 to 3 and 

h±{a) = l±[^^) (5.5) 

For convenience we choose to rewrite the metric of the transverse coordinates as 

dcj2 + a'dnl = dp'' + dy' + p'dnj^^^ + y'dOl^^, (5.6) 

The Hawking temperature in this coordinate system is given by 

T=~ (5.7) 

Assuming that the position of the brane in the y direction is a function of the radial 
coordinate only, the induced metric on the brane is given by 

+ [^)\{i + y\p?)dp' + p'dni_,) (5.8) 

where i and j run from 1 to 3 — A; and a' = p' + y'. One may now write the DBI action 
for the D(7 — 2k) brane and solve the resulting equation of motion to get the profile of the 
brane, from which we can compute the free energy in the field theory. We note that in the 
case of field theories dual to D3 brane geometries the resulting boundary theory is always 
conformal so the equation of motion for the dilaton is solved by a constant for all of the 
defect theories we consider. 

5.1 Brane Embedding 

From the induced metric of a D(7 — 2k) brane (5.8) we can write the DBI action of Nf 
branes including a worldvolume gauge field living on them in the same way as for the 
Sakai-Sugimoto model above 



5D(7-2fc) = rfcAr^/314Vol(Defect) j dpp^-'^h/ J^i^ + v'^) + + 2^a'F)2 (5.9) 
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Note that is the tension of the branes and f3 is the inverse temperature. Vk is the volume 
of an S^~^. Since the action only depends on the derivative of the gauge potential and not 
the gauge potential itself there is an integral of motion 



dC 
dF 



CF 



cf 



27ra'TkNfVk 



CF 



3 — k 

[B + 2'Ka'F) p^-^h^ 



^(l + y'2) + (i5 + 2WFy 



(5.10) 



hence: 



{B + 2TTaFf 



hi cl (1 + y'') 

h+ p2(3-fc)/j3-/c_g2 



(5.11) 



The other equation of motion is 



d_ 
dp 



12 



— Q — — n — 

ap ' op 



dp 



2 G{y',y,p) 



hi 



2y'p^-^h^^'' 

Q {y',y,p) 



(5.12) 



with 



G (y',y,p) 



I hi 



(1 + y'2) + {B + 27ra'FY 



(5.13) 



We find solutions for the embedding function y{p) by using equation (5.11) to eliminate 
(B + 2TTa'F) from (5.12) and solving numerically. For Minkowski embeddings, one can 
show it is not possible to turn on a non constant gauge field A however all chemical 
potentials are accessible via the asymptotic value of a constant gauge field. This means that 
the action, which only depends upon the derivative of A does not depend on the chemical 
potential for Minkowski embeddings, similarly to the Sakai-Sugimoto case. Furthermore, 
for black hole embeddings there is a constraint on cf that comes from (5.11). To ensure 
the reality of {B + 2TTa'F) and hence to ensure the chemical potential is purely imaginary 
we must have 



CF < 2 2 p^ 



(5.14) 



Since we expect there to be Roberge- Weiss lines at ~ exploring the behaviour of the 
theory for small cp is sufficient. 
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Once solutions are found we evaluate the renormalised DBI action'^, dual to the grand 
canonical ensemble in the gauge theory and observe the transition between the two types 
of solution. The chemical potential is the asymptotic value of the worldvolume gauge field 
on the flavour branes, while the expectation value of the phase of the Polyakov loop is 
similarly related to the NS - NS B held. 



a — (Sfij 



Consider the theory when the chemical potential is small. Expanding (5.11) in powers 
oi cp and integrating we find 

a- pjJLi ^^CF + {c%) (5.16) 

where 

^ = /3 / Wi + y'V (5.17) 

J Ph 

^ of course depends on the brane embedding but crucially does not depend on cp. We 
may invert equation (5.16) and substitute into the DBI action (5.9). We find 



-50(7-2^) = rA^//3Vol(53-'=)Vol(Defect) j dp^ 



3-k 



$2 



X ( 1 + ^^2(3-fc)/,3-fc (« - + 0{a- /3/X7)M (5.18) 



This shows that when a — f3pi is small the first order correction to the grand canonical 
ensemble is quadratic in the chemical potential, as expected. 

5.2 The k = 1 Case 

A codimension k = 1 defect means that the glue sector of the theory exists in the full 3 + 1 
dimensions of the boundary while the fundamental fermionic degrees of freedom live on 
a 2 + 1 dimensional subsurface. For a codimension 1 defect in the limit of zero chemical 
potential the phase transition between Minkowski and black hole embeddings is of first 
order as shown in figure 9. One can immediately see the characteristic swallowtail shape 



"The counterterm that exists for all defect theories is Scti = -riV//?Vol(5''^"'")Vol(Defect)^|3^. In the 
case of fe = 2 defects there is an additional counterterm we will describe below 
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Figure 9: For k = 1: Left; The action as a function of ^ for various cf, for black hole embeddings 
only. Right; The action as a function of ^ for dp — 0, both Minkowski and black hole embeddings. 

of the first order phase transition in the grand canonical ensemble for all values of cp and 
therefore for all chemical potentials. In [9] it was shown that there is a first order phase 
transition below some critical value of the chemical potential, above which the transition 
becomes of the second order. In the imaginary chemical potential case the phase transition 
is always of the first order because the Roberge- Weiss lines appear at order so only the 
small imaginary chemical potential behaviour of the phase diagram is accessible. 

One can compute the phase diagram of the theory by examining the point at which the 
single line from the Minkowski embeddings intersects the lines that arise from the black 
hole embeddings. For non zero dp there is no longer a Minkowski embedding allowed, so 
the point where the black hole embeddings for some cp intersect the Minkowski embedding 
line for cp = is taken to be the position of the phase transition. Recall that one may 
turn on an arbitrary chemical potential in the Minkowski embeddings by turning on a 
constant time component of the worldvolume gauge field on flavour branes which does not 
contribute to the DBI action. Therefore, Non zero chemical potentials are accessible in the 
Minkowski phase, while non zero cp is not. 

When we compute the phase diagram we find figure 10. For small fij we the phase 
transition line is given by 

— = 0.298 + 0.273 {a - [Smf + O {a - [Smf (5.19) 
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Figure 10: The first order meson melting transition line in a codimension 1 defect tlieory with 
imaginary chemical potential. 

5.3 The k = 2 Case 

In this setup the fundamental fermions live in a 1 + 1 dimensional subspace of the boundary. 
The gravity theory is quite different in this case because the embedding function diverges 
logarithmically at the boundary, which implies by (5.11) that the worldvolume gauge field 
diverges logarithmically at the boundary as well. We still interpret the constant value of 
the asymptotic embedding function as the quark mass and constant value of the asymptotic 
worldvolume gauge field as the chemical potential. This means the leading order contri- 
bution both of our bulk fields at the boundary is no longer the quark mass and chemical 
potential in the case of the embedding function and worldvolume gauge field respectively, 
but the logarithmic term, the coefficient of which is related to the vacuum expectation 
value of the quark bilinear, and the integral of motion that was defined in (5.10). 

lim y[p) m + clog p + . . . lim At{p) ^ fi + cp log p + . . . (5.20) 

p— ^-oo p— ^oo 

One may show that, as a result of the logarithmic divergences of the embedding function 
and the worldvolume gauge field, the action requires an additional logarithmic counterterm 
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to renormalise it correctly. The additional counterterm is 

Sct2 = -]^TNfP\o\{S'^-'')\o\{Deieci) (c^ + cl) log A 



(5.21) 



When we compute the renormalised action as a function of ^ we find a very similar graph 
to those found before (see figure 11). The discontinuity in the grand canonical ensemble 




Figure 11: The DBI action as a function of ^-^ . On the left is the Minkowski (red) and black 
hole (black) embedding for cp = 0. On the right is a zoomed in version also including black hole 
embeddings with non zero cp (green and blue lines). 

is clearly first order, and when we compute the position of the first order phase transition 
we find figure 12. 

T 



m 



0.2997 + 0.0037 (q - pinf + O {a - I3iiif 



(5.22) 



The phase diagram for codimension 2 systems in ^u^ is still not fully known. It was found 
in [9] that there is a second order phase transition meeting the deconfinement transition so 
there must be a point at some positive /i^ where the phase transition becomes first order, 
since we have shown the phase diagram to be analytic in ^j? for /x^ close to zero. 

5.4 Pressure due to Flavours for Zero Quark Mass 

Similarly to the Sakai-Sugimoto case, the pressure is computed by integrating over the bulk 
directions not on the boundary. 

P = - J (ipd$l3_fc£DBI 



(1 + y'2) + [B + 27ra'Fy 



(5.23) 
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Figure 12: The meson melting pliase transition line of fc = 2 defect theory in the (/i/, — ) plane. 



Dimensional analysis tells us that [P] = 4 — k which is expected because the flavours exist 
only on the defect so the pressure should have the dimension of — [/3 Vol (Defect)] = 4 — k. 
Taking the quarks to be massless simplifies things considerably, because this condition 
implies y' = y = 



NfTkVk I dpp^-^ I 1 + 



3-k 

h 



Using the integral of motion (5.10) 



+ {B + 2TTa'Fy 



(5.24) 



P = -NfTkVk / dpp^-'' ( 1 + ^ 



ph 



3-k 

4\ 



4'" 



(5.25) 



We cannot do this integral analytically, so we expand in powers of the density cp which, 
when we restore the correct factors of the brane tension and a' has dimension 3 — A:, as we 
expect. The leading contribution to the pressure is 

NfTkVkp' 



dp- 



Ph 



[p' + O'^ 



(5.26) 



Which, after integration and renormalisation is 

P 



21+"^ NfTkVkaf-' 



A-k 



(5.27) 
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We can translate this into field theory quantities using (5.7) and the following relations 

1 i?^ 

The volume of a 3 — /c sphere is given by 

therefore the leading contribution to the pressure is given by 

P = 5-^- (5.30) 

(4 - A;)r (1 + 3^) ^ ^ 

The extra symmetry in this model means the leading order expression for the pressure must 

be some dimensionless function times T^~^^ since the temperature is the only dimensionful 

scale in the problem when the chemical potential is zero. The first order correction when 

the chemical potential is non zero is 

^ ^ (5.31) 



{27ra>f {k - 2)NfnVk 



Replacing with field theory quantities as before, the pressure is 

{k - 2)Nf "^^^ ^^■'^^> 

The factor of /c — 2 in the denominator is potentially worrisome, however we shall see that 
when we rewrite the correction in terms of the chemical potential rather than the density 
this factor cancels and the expression is well behaved for all values of k we consider. Using 
the expansion (5.16), with the additional condition that the quarks are massless and we 
find 

\2 2) 

Therefore the first order correction to the pressure is 

4'=-3(/fc-2)A^/7r^(*^-3)Ai 



T^-V^ (5.34) 



r(i-l) 

and the pressure for small density is 



22'=-7iV/7ri('=-3);^2-| / ^ 2(/fc - 2) 



Note that when k = 2 the first order correction to the pressure vanishes. 
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6. Discussion 

We have computed the phase diagrams of the Sakai-Sugimoto model and codimension k 
defect theories for k = 1 and k = 2 and we have shown that these models are analytic in 
/J? when fi'^ is small as is to be hoped since we expect these models to be invariant under 
charge conjugation symmetry, or equivalently the transformation fi — )■ —fi. 

In the models we have considered in this paper we find that in the {fij , T) phase 
diagram all transition lines are first order, meaning that the points where the Roberge- 
Weiss lines meet the phase transition line between low and high temperature phases are 
triple points. The temperature at which the meson melting transition occurs is proportional 
to the square of the imaginary chemical potential, which is similar to what was found in 
[1] for the D3/D7 system. 

One may compare the phase diagrams of codimension k defect theories and we note 
the following, the phase diagram is always of the form 

- = To{k)+T2{k){a-l3fiif + ... (6.1) 
m 

where the coefficients depend on the detail of the defect. We have shown that To is always 
of the same order, but T2 changes significantly between models. Indeed, we note that 
T2 {k + 1) is approximately 100 times smaller than T2 (k), which means the meson melting 
transition line is flatter the higher the codimension of the defect. 

When we compute the pressure for small density the defect theories yield what would 
be expected from dimensional analysis and symmetry arguments - pressure is proportional 
to T^~^ multiplied by a function of The Sakai-Sugimoto model on the other hand does 
not behave as expected, since the pressure goes like at leading order. 
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